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Breathing Modes in Rotating Bose-Condensed Gas: An Exact Diagonalization Study 
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We present an exact diagonalization study of the breathing mode collective excitations for a 
rotating Bose-Einstein condensate of N = 10 spinless bosons interacting via repulsive finite-range 
Gaussian potential and harmonically confined in quasi-two-dimension. The yrast state and the 
low-lying excited states are variationally obtained in given subspaces of the quantized total angular 
momentum L employing the beyond lowest Landau level approximation in slowly rotating regime 
with 0 < L < 2 N. For a given L, the low-energy eigenspectra (bands) are obtained in weakly to 
moderately interacting regime. Further, for a given interaction, the split in low-lying eigenenergies 
with increasing L is the precursor to spontaneous symmetry breaking of the axisymmetry associated 
with the entry of the first vortex. With increase in repulsive interaction, the value of the first 
breathing mode increases for stable total angular momentum states L = 0 and N, but decreases for 
intermediate 0 < L < N metastable states. The position of the observed first breathing modes in 
the eigenspectrum remains unchanged as the interaction is varied over several orders of magnitude. 
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I. INTRODUCTION 


II. THEORETICAL FORMALISM 


Ever since the experimental realization of Bose- 
Einstein condensation (BEC) in ultra-cold alkali atomic 
vapours in a harmonic trap [Ml, the study of collective 
excitations in such systems has been an important sub¬ 
ject of research in quantum many-body physics H M- 
The decisive experimental control over the density, the 
effective dimensionality and the atom-atom interaction 
strength [Bj, makes these systems an outstanding one to 
study subtle quantum many-body effects, in particular, 
collective excitations. 

The breathing mode (monopole oscillation or uniform 
radial expansion and contraction) is one of the most im¬ 
portant collective excitations Eli used as a diagnostic 
tool for many-body effects [9j. In recent years, a num¬ 
ber of theoretical studies have demonstrated that the 
quantum breathing mode is ideally suited to estimate 
the atom-atom coupling strength, its kinetic and interac¬ 
tion energies and other such observables [lOj, in a trapped 
atomic vapour system. This leads to a novel kind of spec¬ 
troscopy of trapped systems. 

For few-body systems realized in lower dimension 
[nun!, the behaviour of breathing modes can be stud¬ 
ied to a very high degree of precision mm ■ An under¬ 
standing of the physics of few-body systems may then be 
extrapolated to larger systems to gain an insight into the 
beyond mean-field physics of macroscopic ensembles. In 
this note, we follow this approach and present an exact 
diagonalization study of the many-body effects of inter¬ 
action and rotation (with quantized L) on the dynamics 
of breathing modes. 
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We consider a system of N interacting spinless bosons 
each of mass M , trapped in a harmonic potential V (r) = 
\M + u> 2 z 2 ). The trap is subjected to an external 

rotation about the z-axis with angular velocity fl = f le z . 
The trapping potential V (r) is assumed to be highly 
anisotropic with A 2 = w 2 /w_l 1, so that the many- 
body dynamics along z-axis is frozen. The system is thus 
effectively quasi-two-dimensional (quasi-2D) with x-y ro¬ 
tational symmetry. Choosing hui± and aj_ = yjfi/M uj±_ 
as units of energy and length respectively, the many-body 
Hamiltonian in the co-rotating frame is given as 


N r 


j =1 




1 N 


(!) 

where f l = fl/ui± (< 1) is the dimensionless angular ve¬ 
locity and L (scaled by K) is the z projection of the total 
angular momentum operator. The inter-particle interac¬ 
tion U is described by the Gaussian potential 


U (I r i-rj|) = 


g 2 

27T(J^ 


exp 


(r±i - r±j) 

2o\ 


(5 (zi - Zj) 


( 2 ) 

with width a± (scaled by aj_) being the effective range 
of two-body interaction. The dimensionless parameter 
g 2 = 47ra s /a_L measures the strength of the two-body in¬ 
teraction with a s being the s-wave scattering length for 
particle-particle collision. We assume that the scatter¬ 
ing length is positive (a s > 0) so that the effective finite- 
range interaction is repulsive. The above finite-range 
Gaussian interaction potential is expandable within a fi¬ 
nite number of single-particle basis and hence computa¬ 
tionally feasible [HI flBl . In the limit a± —> 0, the Gaus¬ 
sian potential in Eq. (pj) reduces to the usual zero-range 
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contact potential g 2 S (r t — r j). 

To obtain the eigenenergies and the corresponding 
eigenstates, we employ exact diagonalization of the 
Hamiltonian matrix in different subspaces of L using 
Davidson algorithm na with inclusion of higher Landau 
levels in constructing the many-body basis states Tl5] . 
It is to be noted that for a many-body system under 
consideration here, the characteristic energy scale for the 
interaction is determined by the dimensionless parameter 
Na s /a±. Owing to the increasing dimensionality of the 
Hilbert space with N making the computation imprac¬ 
tical, we vary a s so as to achieve the value of Na s /a± 
relevant to experimental situation [5]. 


III. RESULTS AND DISCUSSIONS 


We consider a system of N = 10 bosonic atoms of 8l Rb 
in a quasi-2D harmonic trap with confining frequency 
w_l = 2-7T x 220 Hz and the z-asymmetry parameter 
A 2 = u) z /w_l = y/8. The condensate has extension 
a z = Ti/Mui z in the z-direction and its dynamics 
along this axis is taken to be completely frozen. Recent 
advancements in atomic physics has made it possible 
to tune the low-energy atom-atom scattering length in 
ultra-cold atomic vapours using Feshbach resonance [B]. 
Accordingly in the calculations presented here, the pa¬ 
rameters of the two-body interaction potential in Eq. (|2j) 
has been chosen as a± =0.1 and a s = 10 ao, 100 ao and 
1000 do where, ao = 0.05292 nrn is the Bohr radius. The 
corresponding dimensionless parameter g 2 = 47ra s /a_L 
defined above turns out to be 0.009151, 0.09151 and 
0.9151, respectively. 

Following m, we examine here the breathing modes 
in a rotating system of N = 10 bosons within the 
quantized total angular momentum regime 0 < L < 2TV. 
In Fig. [I] we present the low-energy eigenspectra 
(bands) for different total angular momentum L states 
with three different values of interaction parameter g 2 
(parametrized by the scattering length a s ). Eigenstates 
having the same total angular momenta constitute a 
L series (or bands). The ith eigenstate of the L series 
is denoted as Li and the corresponding eigenenergy 
as E{Li). We observe that corresponding to each L, 
the quasi-degenerate low-lying eigene nergie s spiff to 
form energy bands for weakly, Figs. l(a)| and |l(b)| 
to moderately interacting regime, Fig. l(c)| and the 
respective energy gaps associated with the low-lying 
eigenstates are reduced. We also observe that for a 
given interaction strength, the low-lying eigenenergies 
further split with increase in L to fill the energy gaps, 
as seen in Figs. |l(b)| and 1(c) This split in low-lying 
eigenenergies with increasing L leads to spontaneous 
symmetry breaking of the axisymmetry associated with 
the entrance of the first vortex [T9] in the angular 
momentum regime 0 < L < N. 

To analyze the breathing modes, we present in 
Table [T] the low-lying eigenenergies E {Li) corresponding 
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(b) a s = 100 ao 



FIG. 1. For N = 10 condensed bosons in a quasi-2D harmonic 
trap, the low-lying energy eigenspectrum (in units of hu>±) 
versus the total angular momentum L (scaled by h) for three 
different values of repulsive interaction g 2 (parametrized by 
a s ) and a± = 0.1 in Eq. In the present calculation, up to 
thirty low-lying energy eigenstates have been found. 


to angular momentum states 0 < L < 11 with the 
interaction parameter value g 2 = 0.09151. The lowest 
energy state L\ corresponding to an L is referred to 
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TABLE I. (Color online) The eigenenergies of the L; states for N = 10 bosons with interaction parameters g 2 = 0.09151 and 
a± = 0.1 in the total angular momentum regime 0 < L < 11. The L\ states are the yrast states (or ground modes) and the L; 
states O2, I2, 23, 34, 45, 56, 6s, 79 , 8n, 9i3, IO24, H27 are the observed first breathing modes. 


i 

0 

11 

L = 1 

L = 2 

CO 

II 

1-4 

L = 4 

L = 5 

L = 6 

II 

00 

II 

L = 9 

L = 10 

L = 11 

1 

24.5663 

25.5663 

26.5217 

27.4976 

28.4778 

29.4524 

30.4299 

31.4062 

32.3822 

33.3583 

34.3339 

35.3340 

2 

26.5664 

27.5217 

26.5665 

27.5216 

28.4963 

29.4774 

30.4416 

31.4309 

32.4039 

33.3830 

34.3578 

35.3459 

3 

26.5688 

27.5447 

28.4996 

27.5670 

28.5202 

29.4952 

30.4545 

31.4361 

32.4102 

33.3864 

34.3630 

35.3563 

4 

28.4833 

27.5668 

28.5218 

29.4781 

28.5675 

29.5184 

30.4752 

31.4545 

32.4308 

33.4009 

34.3746 

35.3582 

5 

28.5218 

27.5687 

28.5231 

29.4977 

30.4540 

29.5679 

30.4940 

31.4716 

32.4321 

33.4109 

34.3819 

35.3648 

6 

28.5224 

29.4607 

28.5392 

29.4990 

30.4759 

31.4315 

30.5164 

31.4926 

32.4522 

33.4261 

34.3830 

35.3721 

7 

28.5438 

29.4832 

28.5671 

29.5051 

30.4779 

31.4440 

30.5681 

31.5143 

32.4672 

33.4313 

34.3846 

35.3769 

8 

28.5448 

29.5007 

28.5680 

29.5204 

30.4807 

31.4523 

32.4074 

31.5678 

32.4911 

33.4482 

34.3948 

35.3819 

9 

28.5662 

29.5018 

30.4454 

29.5228 

30.4953 

31.4557 

32.4285 

33.3834 

32.5119 

33.4627 

34.4009 

35.3826 

10 

28.5676 

29.5221 

30.4598 

29.5349 

30.4982 

31.4619 

32.4324 

33.4047 

32.5671 

33.4894 

34.4056 

35.3848 

11 

28.5687 

29.5228 

30.4787 

29.5661 

30.5002 

31.4734 

32.4357 

33.4060 

34.3591 

33.5091 

34.4071 

35.3921 

12 

30.4398 

29.5233 

30.4798 

29.5683 

30.5184 

31.4759 

32.4399 

33.4106 

34.3808 

33.5658 

34.4144 

35.3986 

13 

30.4617 

29.5365 

30.4819 

31.4223 

30.5205 

31.4779 

32.4454 

33.4138 

34.3848 

35.3347 

34.4225 

35.4007 

14 

30.4627 

29.5416 

30.4990 

31.4448 

30.5312 

31.4903 

32.4529 

33.4257 

34.3876 

35.3601 

34.4279 

35.4064 

15 

30.4833 

29.5444 

30.4998 

31.4557 

30.5647 

31.4960 

32.4549 

33.4306 

34.3915 

35.3635 

34.4308 

35.4127 

16 

30.4836 

29.5659 

30.5006 

31.4578 

30.5686 

31.4977 

32.4579 

33.4322 

34.4033 

35.3718 

34.4344 

35.4203 

17 

30.4944 

29.5676 

30.5051 

31.4591 

32.4027 

31.5161 

32.4680 

33.4331 

34.4051 

35.3814 

34.4470 

35.4268 

18 

30.4970 

29.5684 

30.5153 

31.4771 

32.4137 

31.5171 

32.4724 

33.4371 

34.4097 

35.3831 

34.4530 

35.4304 

19 

30.5037 

31.4247 

30.5209 

31.4779 

32.4251 

31.5278 

32.4753 

33.4407 

34.4113 

35.3867 

34.4723 

35.4317 

20 

30.5040 

31.4407 

30.5224 

31.4791 

32.4330 

31.5629 

32.4849 

33.4480 

34.4127 

35.3883 

34.4858 

35.4450 

21 

30.5212 

31.4456 

30.5231 

31.4803 

32.4353 

31.5687 

32.4932 

33.4530 

34.4159 

35.3914 

34.4936 

35.4517 

22 

30.5219 

31.4601 

30.5238 

31.4815 

32.4429 

33.3831 

32.4956 

33.4552 

34.4275 

35.3954 

34.5161 

35.4699 

23 

30.5228 

31.4614 

30.5339 

31.4916 

32.4468 

33.4038 

32.5127 

33.4624 

34.4296 

35.3993 

34.5672 

35.4854 

24 

30.5229 

31.4625 

30.5383 

31.4957 

32.4534 

33.4078 

32.5141 

33.4677 

34.4305 

35.4033 

36.3342 

35.4926 

25 

30.5242 

31.4676 

30.5442 

31.4980 

32.4566 

33.4090 

32.5245 

33.4709 

34.4330 

35.4069 

36.3346 

35.5146 

26 

30.5354 

31.4769 

30.5642 

31.4986 

32.4571 

33.4129 

32.5608 

33.4794 

34.4367 

35.4092 

36.3466 

35.5669 

27 

30.5373 

31.4798 

30.5662 

31.5006 

32.4585 

33.4253 

32.5683 

33.4899 

34.4407 

35.4098 

36.3558 

37.3287 

28 

30.5416 

31.4822 

30.5689 

31.5040 

32.4619 

33.4262 

34.3632 

33.4937 

34.4490 

35.4123 

36.3572 

37.3336 

29 

30.5426 

31.4834 

32.4033 

31.5088 

32.4665 

33.4297 

34.3816 

33.5085 

34.4519 

35.4203 

36.3586 

37.3344 

30 

30.5589 

31.4912 

32.4232 

31.5188 

32.4749 

33.4340 

34.3843 

33.5114 

34.4572 

35.4246 

36.3599 

37.3352 

as 

the yrast 

state (we would occasionally call it 

the 

E( 10 24 ) - 

E( 10i) 

= 2.0003 

for L 

= 10. Thus, the 


ground mode). We observe that for the non-rotating 
L = 0 states, E( O 2 ) — -E7(0i) = 2.0001 (in units of hw±). 
The states 0i and 0 2 are, thus, taken to be the ground 
mode (the yrast state) and the first breathing mode, 
respectively, for L = 0. It was pointed out by Pitaevskii 
and Rosch [7 that a purely 2D bosonic system with 
zero-range interaction exhibits breathing modes arising 
from 50(2,1) symmetry. The energy difference between 
adjacent breathing mode levels is 2Tluj±. Although in 
the present work, the finite-range Gaussian interaction 
potential has been used to replace the zero-range 
(d-function) interaction, the feature of 2huj_ spacing is 
seen to persist in the energy eigenspectrum. Further, 
the 2huj± spacing also appears for the rotating L > 0 
states. For example, we find E( 1 2 ) — I?(li) = 1.9554 
for L = 1, E{ 5 6 ) - E( 5i) = 1.9791 for L = 5 and 


breathing modes existing in 2D systems with zero-range 
(5-function) interaction potential may also exist in our 
quasi-2D system with Gaussian interaction potential. 

For the non-rotating L = 0 case (column 2 in the 
Table [I]), the position of the first breathing mode 0 2 is 
always immediately above the yrast state 0i. However 
as L increases, the positions (labelled by the index i 
in column 1) of the first breathing modes are shifted 
to higher eigenstates, forming a stair-like pattern of 
eigenenergies corresponding to the breathing modes. 
We further observe from the Table |T] that the shifts in 
the positions of the first breathing modes are small for 
states with L < N; however, for the first vortex state 
with L = N = 10, we notice an appreciable shift in the 
position of the first breathing mode. It might be due 
to rearrangement of bosons on the nucleation of first 
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central vortex. We have found that as the repulsive 
interaction is increased over three orders of magnitude, 
the positions of the first breathing modes in the energy 
eigenspectra for a given L, remain unchanged. 

We also observe the second and third breathing 
modes with approximately and 6huj± spacing. For 

example, the states O4, lg, 2 g , 3i3, 447, 522 and 62 s are 
the second breathing modes. Similarly, the states O12, 
I 19 and 229 are the third breathing modes. Moreover as 
is the case for the first breathing mode, for a given L, 
the positions of these higher breathing modes are also 
found to be independent of the inter-particle interaction 
strength. 

TABLE II. For N = 10 condensed bosons, the values of 
the breathing modes Ebm (in units of hco±) i.e. the differ¬ 
ence between the eigenenergies corresponding to the breathing 
state (Li) and the yrast state (L 1 ), for three different values 
of the dimensionless interaction parameter g 2 (parametrized 
by a a ) with Gaussian width a± = 0.1 in Eq. d2p. With in¬ 
crease in repulsive interaction, the value of the first breath¬ 
ing mode increases for stable total angular momentum states 
L = 0 and N, but decreases for intermediate 0 < L < N 
metastable states. 


L 

E 

BM 

Ebm (10a o ) E B m (100a 0 ) E) 

bm (1000a 

0 

E( 0 a ) 

-E( 0i) 

1.9999 

2.0001 

2.0111 

1 

E( D) 

~E{ li) 

1.9952 

1.9554 

1.7837 

2 

E( 2 3 ) 

~E( 2 r) 

1.9978 

1.9779 

1.8062 

3 

E( 3 4 ) 

~E( 3i) 

1.9978 

1.9805 

1.9357 

4 

E( 4s) 

-E( 4i) 

1.9976 

1.9762 

1.8287 

5 

E( 5 6 ) 

~E{ 5r) 

1.9978 

1.9791 

1.8542 

6 

E( 6 8 ) 

-£(61) 

1.9976 

1.9775 

1.8620 

7 

E( 7 9 ) 

~E{ 7i) 

1.9976 

1.9772 

1.8403 

8 

E( 8n) 

-£(8r) 

1.9976 

1.9768 

1.8331 

9 

E(9 13) 

-£(9i) 

1.9976 

1.9764 

1.7263 

10 £(10 24 ) 

-£(10i) 

2.0000 

2.0003 

2.0189 

11 £(11 27 ) 

-£(lli) 

1.9994 

1.9947 

1.9792 


The results of the first breathing modes in subspaces of 
quantized L are presented in Table |TT] for three different 
values of the interaction parameter g 2 . We observe 
that for a given L , the values of the first breathing 
mode calculated with weak interactions corresponding 
to the choice g 2 = 0.009151 and 0.09151 are very close 
(column 3 and 4) but significantly different from the ones 
calculated for moderate interaction corresponding to the 
choice g 2 = 0.9151 (column 5). For the non-rotating 


L = 0 and the single-vortex L = N = 10 states, the 
values of the breathing modes increase albeit very 
weakly with increase in repulsive interaction. On the 
contrary for the intermediate 0 < L < N metastable 
states, the values of the breathing modes are appreciably 
different from the non-rotating and single-vortex states 
and decreases significantly with increase in repulsive 
interaction. 

IV. CONCLUSIONS 

In conclusion, we have examined the breathing mode 
collective excitations of rotating Bose-Einstein conden¬ 
sate with finite-range Gaussian interaction in a quasi-2D 
harmonic trap. By exact diagonalization of the many- 
body Hamiltonian matrix in beyond lowest-Landau-level 
approximation, the eigenenergies of the ground and the 
low-lying excited states are obtained with total angular 
momenta 0 < L < 2TV, corresponding to slowly rotating 
regime. We have presented the low-energy eigenspec¬ 
tra (bands) of the rotating system for three representa¬ 
tive values of different interacting regimes in subspaces of 
quantized L. For a given L, the degeneracy of eigenener¬ 
gies is lifted (forming energy bands) for weakly to mod¬ 
erately interacting regime and correspondingly the en¬ 
ergy gaps associated with the low-lying eigenstates are 
reduced. We find that the repulsive interaction influ¬ 
ences the values of the first breathing modes in two dif¬ 
ferent ways. For stable L = 0, N states the value of the 
first breathing mode increases whereas for intermediate 
0 < L < N metastable states it decreases with increase 
in repulsive interaction. For metastable states with an¬ 
gular momentum L < N, we observe a small shift in the 
position of the first breathing mode in the eigenspectrum 
whereas for the vortex state (L = N), the shift becomes 
significant. The position of the observed first breathing 
mode in the eigenspectrum remains unchanged as the in¬ 
teraction is varied over several orders of magnitude. 

ACKNOWLEDGMENTS 

The present work was put forth during the 4th Interna¬ 
tional conference on “Current Developments in Atomic, 
Molecular Optical and Nano Physics with Application” 
CDAMOP-2015, at Department of Physics and Astro¬ 
physics, University of Delhi, India. The authors are 
thankful to the organisers for providing a platform for 
exchange of ideas. 


[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C.E. 
Wieman, and E. A. Cornell, Science 269 , 198 (1995). 

[2] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van 
Druten, D. S. Durfee, D. M. Kurn, and W. Ketterle, 
Phys. Rev. Lett. 75 , 3969 (1995). 


[3] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. 
Hulet, Phys. Rev. Lett. 75 , 1687 (1995). 

[4] D. S. Jin, J. R. Ensher, M. R. Matthews, C. E. Wieman, 
and E. A. Cornell, Phys. Rev. Lett. 77 , 420 (1996). 








5 


[51 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80 , 885 (2008). 

[6] S. Inouye, M. R. Andrews, J. Stenger, H.-J. Miesner, D. 
M. Stamper-Kurn, and W. Ketterle, Nature (London) 
392 , 151 (1998). 

[7] L. P. Pitaevskii, and A. Rosch, Phys. Rev. A 55, R853 
(1997). 

[8] F. Chevy, V. Bretin, P. Rosenbusch, K. W. Madison, and 
J. Dalibard, Phys. Rev. Lett. 88, 250402 (2002). 

[9] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

[10] C. R. McDonald, G. Orlando, J. W. Abraham, D. 
Hochstuhl, M. Bonitz, and T. Brabec, Phys. Rev. Lett. 
Ill, 256801 (2013). 

[11] F. Serwane, G. Zlirn, T. Lompe, T. B. Ottenstein, A. N. 
Wenz, and S. Jochim, Science 332 , 336 (2011). 


[12] E. Haller, R. Hart, M. J. Mark, J. G. Danzl, L. Re- 
ichsollner, M. Gustavsson, M. Dalmonte, G. Pupillo, and 
H.-C. Nagerl, Nature (London) 466, 597 (2010). 

[13] C. G. Bao, Y. Z. He, G. M. Huang, and T. Y. Shi, Phys. 
Rev. A 65 , 022508 (2002). 

[14] M. Olshanii, H. Perrin, and V. Lorent, Phys. Rev. Lett. 
105 , 095302 (2010). 

[15] J. Christensson, C. Forssen, S. Aberg, and S. M. 
Reimann, Phys. Rev. A 79 , 012707 (2009). 

[16] R. A. Doganov, S. Klaiman, O. E. Alon, A. I. Streltsov, 
and L. S. Cederbaum, Phys. Rev. A 87, 033631 (2013). 

[17] E. R. Davidson, J. Comput. Phys. 17, 87 (1975). 

[18] M. A. H. Ahsan, and N. Kumar, Phys. Rev. A 64, 013608 

( 2001 ). 

[19] T. Nakajima, and M. Ueda, Phys. Rev. Lett. 91, 140401 
(2003). 



